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Coherent Anti-Stokes Raman Scattering (CARS) processes are "coherent," but the phase of the 
anti-Stokes radiation is usually lost by most incoherent spectroscopic CARS measurements. We 
propose a novel Raman microscopy imaging method called Nonlinear Interferometric Vibrational 
Imaging, which measures Raman spectra by obtaining the temporal anti-Stokes signal through 
nonlinear interferometry. With a more complete knowledge of the anti-Stokes signal, we show 
through simulations that a high-resolution Raman spectrum can be obtained of a molecule in a 
single pulse using broadband radiation. This could be useful for identifying the three-dimensional 
spatial distribution of molecular species in tissue. 



Functional imaging techniques have been developed to 
provide insight into biological processes. Optical func- 
tional imaging is frequently limited because dyes or mark- 
ers must be introduced that alter or damage biologi- 
cal tissues. Because it is preferable to use endogenous 
properties of tissues to identify molecular contents, meth- 
ods such as infrared and Raman spectroscopy are used. 
In particular, Coherent Anti-Stokes Raman Scattering 
(CARS) processes have been successfully integrated with 
confocal scanning microscopes Jj, to achieve three- 
dimensional molecular images of vibrational resonances. 
However, existing instruments measure only the total 
power of the received anti-Stokes radiation. We propose 
a novel method. Nonlinear Interferometric Vibrational 
Imaging (NIVI), which utilizes nonlinear interferometry 
to measure the amplitude and phase of the anti-Stokes 
light. An experimental demonstration of the principle 
of this technique has been demonstrated. This addi- 
tional phase information facilitates the inference of the 
amplitude and phase of the nonlinear susceptibility of the 
molecule. By utilizing NIVI with properly designed illu- 
minating radiation, a large region of the amplitude and 
phase of the Raman spectrum can be sampled in a single 
brief pulse. 

Coherent Anti-Stokes Raman Scattering processes 
have only recently been used to probe biological spec- 
imens. The appeal of CARS is that it can probe the 
density of molecules with a particular Raman resonance 
frequency while exposing the specimen to relatively low 
levels of illumination. A typical CARS process illumi- 
nates the specimen with a pump pulse of frequency oji , 
and a Stokes pulse of frequency L02, which are separated 



by the vibrational frequency of interest Q ~ loi — lu2. If 
a molecule with a Raman resonance at frequency Q is 
present, an anti-Stokes pulse at frequency 2lui — L02 is 
produced. In CARS microscopy, tightly focused pump 
and Stokes beams are scanned through the specimen, 
and the anti-Stokes photon count is measured at each 
point. This photon count is proportional to the square 
of the molecular bond density and the magnitude of the 
Raman susceptibility. Nonlinear interferometry has been 
used to characterize the magnitude of the stimulated Ra- 
man scatteringnonlinearity in liquids 5] , and also of two 
CARS signals p in gases. 

The model of CARS with broadband pulses used here 
is similar to that described by 0, ■ We do not assume 
that the illuminating radiation is narrowband. However, 
we stipulate that the molecule does not resonantly inter- 
act directly with any of the frequencies inside the illu- 
mination bandwidth or the generated anti-Stokes band- 
width of the optical signal. The CARS process is com- 
posed of two stimulated Raman scattering processes in- 
volving four photons that begins and ends with the 
molecule in the ground state. To describe CARS, we de- 
note the electric field incident on the molecule as EiiiS). 
The first process is modeled by Eq. ^ and excites the 
nonlinear dipole polarization of the resonant transition. 



p(3)(17) =x^^'^{^) / Ei{u ^^)Ei{uY duj (step 1) (1) 



Eo{uj) = J E^{uJ - n)p^^\n) dn (step 2) (2) 





9. 



Each pair of frequencies that are separated by a reso- 
nance of the molecule at frequency produces a nonlin- 
ear polarization in the molecule. Another way to look 
at Eq. n is that, in the time domain, the molecule has 
a nonlinear polarization that is driven not by the elec- 
tric field but by the instantaneous intensity envelope of 
the signal. In this formulation, we are neglecting any 
changes in x^^^ that are dependent on the "carrier" en- 
velope frequency on which the beats are imposed. There- 
fore any pulse train with intensity beats at the resonance 
frequency will stimulate the nonlinear polarization. Ex- 
amples of this are pulses that are modulated periodically 
in the spectral domain with period B 1 ^^'^ i'^" 
terfering two relatively delayed chirped pulses to achieve 
a beat frequency proprortional to time delay The 
second step creates anti-Stokes radiation by mixing the 
incoming radiation field with the polarization in the time 
domain, and is modeled by Eq.|2| The Eqs.^ and fallow 
one to calculate the emitted CARS radiation Eo{uj) for 
a given Ei{uj) and x^^H^)- While these relations do not 
constitute a Hnear relationship between Ei{(jj) and Eo{uj), 
there is a linear dependence of Eo{^) on x^^^ i^) given an 
input field Ei{uj). This suggests that with knowledge of 
the complex Eo{lu), one can do linear estimation to find 
x('^)(rj). The advantage of NIVI over incoherent detec- 
tion is that nonlinear interferometry enables the recovery 
of the complex Eo{oj)- With a properly designed input 
pulse Ei{uj), the nonlinear susceptibility can be found in 
a particular frequency range. 

NIVI takes advantage of the coherent nature of the 
CARS process to allow the phase-sensitive measurement 
of the anti-Stokes radiation. Conventional linear inter- 
ferometry involves splitting a source light beam into two 
parts, each of which scatters linearly in the field, and 
which are then recombined and detected. NIVI differs in 
that a CARS process occurs to one of the split beams. 

NIVI can be implemented with the setup detailed 
in Fig. n We start with a broadband, phase-locked 
source of light pulses, such as those from a mode-locked 
laser. Sources that can produce such light are ultra- 
broadband mode-locked Ti-sapphire oscillators ^"^^ 
supercontinuum sources 0, • Because the source 
is phase-locked, there is a deterministic relationship be- 
tween phases at various frequencies. This deterministic 
relationship will be preserved by coherent processes such 
as CARS. To utilize this determinacy, the source band- 
width is split into higher and lower frequency bands, as 
shown in Fig. with a dichroic beamsplitter. The higher 
frequency band will be a reference signal R{uj), and will 
correspond to the bandwidth of the anti-Stokes frequen- 
cies produced by a sample. The lower frequency band is 
temporally shaped to stimulate CARS in the sample, a 
signal we denote by the frequency spectrum Eiiuj). Some 
of the illumination signal will be converted to anti-Stokes 
radiation by the sample. Because CARS processes are 
usually rather weak, we will assume that any new anti- 



Broadband 

Continuum 
Pulsed 
Source 



Dichroic 
Low Pass 

Beam 

Splitter 

\ Pump/Stokes 
^^Frequencies 




lempora' 
Pulse 
Shaper 




CARS 




Sample 



High Pass ^ 
Filter 



Reference Frequencies 
(same as anti-Stokes 
bandwidth) 



\ 



Anti-Stokes 
Frequencies Only 

I. 

Power 
Detector 



I. 

Power 
Detector 



50:50 
Beam 
Splitter 



FIG. 1: Schematic of NIVI with a broadband continuum 
source, using the continuum for the reference pulse. The 
detection scheme is simplified and in practice a single-shot 
detection scheme such as TJl could be used. 



Stokes radiation created in the same bandwidth as the 
illumination will be inseparable from the illumination. 
Therefore we will discard all anti-Stokes light inside the 
illumination bandwidth with a high-pass frequency fil- 
ter. The remaining anti-Stokes light that passes through 
the filter, which we denote by the spectrum i?o(w), corre- 
sponds to frequencies in the reference signal. We combine 
the reference signal R{lo) and the anti-Stokes spectrum 
Eo{w) with a 50:50 beam splitter, and utilize balanced de- 
tection to measure the intereference component on two 
photodetectors. There is a delay of time At placed in 
the reference path to facilitate measuring the temporal 
cross-correlation between the reference and anti-Stokes 
signals. The difference in the two intensities A/(Ai) as 
a function of delay between the signal and reference will 
be: 



A/(At) = I+-I- 



4 Re 



|£'o(w)*^(w)exp(iwAi)| duj 



(3) 

If we call I{uj) the Fourier transform of A/(At) with re- 
spect to At, we find that A/(w) = i Eoi^^)* Rioj) ■ Thus 
the measured data retains its linear relationship with re- 
spect to the anti-Stokes spectrum Eo{'-^) and therefore 
the nonlinear susceptibility x^^H^)- 

Besides the ability to find the complex- valued x'^' 
interferometry eliminates the need for photon-counting 
detectors. Another advantage is that interference will 
only occur when the anti-Stokes light and the reference 
light arrive at the beam splitter at the same time. Be- 
cause of this, temporal gating can be used to produce 
three-dimensional vibrational images in a manner analo- 
gous to Optical Coherence Tomography 0, 0, 0| . Co- 
herent detection is also far less sensitive to stray light 
than photon counting. Because of this, NIVI may be 
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more adaptable to various scanning configurations and 
environments outside the laboratory. 

To show that NIVI can measure intervals of the Ra- 
man spectrum in a single pulse, a pulse must be designed 
that can stimulate molecules in a broad Raman spec- 
trum. The approach we take creates beats that instead 
of being of a constant frequency 1^ , will be themselves 



chirped. This can be accomplished by combining two 
chirped pulses with a relative delay, but with different 
chirp rates. If we have a transform-limited pulsed source 
of center frequency loq and bandwidth Auj, and we wish 
to sweep the beat frequency from fi^ to Qh in time T, 
we can design a pulse Ei{oj) such that: 



AH = i?oCos(l(^) 

( 1-) exp 
Ei{uj) — otherwise 
where a = 



(iii^) exp 



-i(uj—ujQ)T 
2 



i(w-wo)' 



for Uq 



2{a+l3) 

At 
2 



(4) 



2T 



and T 



2 \ Aoj-Of: 



The variable a is the common chirp to both pulses, /3 
is the difference chirp, r is the time delay between the 
two pulses, and k is the difference in field magnitude be- 
tween the two pulses. The pulse bandwidth has been 
apodized with a cosine window because in practice it 
seems to help the stability of the inversion. Note that 
the bandwidth of the source Alu must exceed fin so that 
beats can be formed at all Raman frequencies. When 
creating the pulse, the chirp time T will control the reso- 
lution with which one will be able to resolve frequencies 
in the Raman spectrum. The largest practical T is de- 
termined by the dephasing time of the resonances, which 
in most liquids is on the order of picoseconds. 

To demonstrate the feasibility of NIVI, we simulate 
the illumination of a target molecule with the broadband 
pulse of Eq.0]and use the returned signal to estimate the 
complex susceptibility x'^^H^)- We will show two simu- 
lations: one that is able to probe a wide bandwidth of 
Raman resonances in a single pulse, and the other which 
is able to distinguish between two nearby resonances. We 
take as our hypothetical laser source a mode-locked Ti- 
sapphire laser that can produce a pulse with a uniform 
bandwidth from 700-1000 nm, and the setup of Fig. ^ 
The bandwidth from 800-1000 nm will be reserved for 
stimulating CARS, with the remainder used as the ref- 
erence signal. For the first simulation, the CARS excita- 
tion bandwidth will be shaped such that fi^ = 700 cm^^, 
^Ih — 1300 cm^^, and T = 5 ps. To show that the sys- 
tem can reconstruct several simultaneous resonances over 
the entire bandwidth, we create a hypothetical x^^H^) 
with several Lorentzian resonances centered at 800 cm~^, 
900 cm^^, 1000 cm~^, and 1100 cm~^. These frequen- 
cies are in the Raman "fingerprint" region and would be 
useful for practical molecular identification. 

The simulation was implemented by sampling the spec- 
tra of x''^\^), E,{uj), and with 20,000 points 
spaced at equal intervals from cm~^ to 20000 cm^^ in 



1.0 cm~^ steps. The cross-correlations of Eqs. ^ and [21 
were computed using the Fast Fourier Transform. These 
two equations form a "forward" CARS linear operator 
computing x^^'^ {^) from Eo{ijj) which we call A.{ijj^ fi). To 
find the inverse of this operator, we used the Tikhonov- 
regularized least-squares inversion operator, which is for- 
mally denoted by A* = (At A + el)-i A^. The Tikhonov 
regularization was included to improve the stability of the 
inverse and to account for potential noise sources such 
as thermal and photon noise in a practical experiment. 
The constant e is chosen to account for the magnitude of 
additive white Gaussian noise in a realistic experiment. 
In practice, A* was computed using the preconditioned 
conjugate gradient |l8j method to avoid the very com- 
putationally expensive direct matrix inversion. While we 
do not model a real noise source here, Tikhonov regular- 
ization adjusts the inverse operator such that features of 
the estimated spectrum x^^H^) that would be unstable 
due to insufficient information for reconstruction would 
tend towards zero. 

The left column of Figure |21 shows the temporal and 
spectral shapes of the input pulse. Part (a) shows two 
chirped pulses that partially overlap, producing a beat 
pattern that stimulates the resonance. Part (b) shows 
the power spectrum of the input pulse. Part (c) shows 
the anti-Stokes radiation spectrum, that is calculated us- 
ing Eqs.HandlU Because the excitation light is assumed 
to be much more powerful that the anti-Stokes light, we 
filter out all of the excitation bandwidth and utilize only 
wavelengths shorter than 800 nm for the inverse. The 
right column of Figure [21 shows the original and recon- 
structed x^^^ iP) ■ Pa-rt (d) is the magnitude of the spec- 
trum of the intensity of the original pulse, i.e. the beat 
frequencies of the pulse. It shows the possible measurable 
Raman frequencies with this pulse. Part (e) shows the 
original x^^^ (^) Raman spectrum magnitude. Finally, 
part (f) is the Tihkonov- regularized least-squares recon- 
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FIG. 2: NIVI input and output pulses, and the original and 
reconstructed Raman spectra in first simulation, (a) Tem- 
poral shape of amplitude of input pulse, (b) Power spectral 
density of input pulse, (c) Power spectral density of out- 
put pulse, (d) Power spectrum of beat frequencies of input 
pulse, (e) Magnitude of Raman susceptibility of hypothet- 
ical molecule, (f) Magnitude of least-squares reconstructed 
Raman susceptibility of hypothetical molecule. 



structed x^^^(f2) based on only the anti-Stokes frequen- 
cies from 700-800 nm. In simulation, all of the spectral 
lines can be recovered. The minimum discernible sep- 
aration in Raman frequencies tends to increase as the 
Raman frequency decreases because the anti-Stokes ra- 
diation created by lower frequency resonances tends to 
overlap the original spectrum more. 

As a second demonstration with two closely spaced 
Raman lines, we consider deoxyribonucleic acid (DNA), 
which would be contained in the nucleus of a cell, and ri- 
bonucleic acid (RNA) located throughout the cell. Both 
macromolecules have PO2 phosphodiester resonances, 
but the resonance occurs in DNA at 1094 cm~^ and in 
RNA at 1101 cm"-'^. To show that a properly designed 
pulse can recover both resonances distinctly, we design 
a x'^'(^) that has resonances both at 1094 cm~^ and 
1101 cm~^, that could be created by mixing DNA and 
RNA. To probe this mixture, we create a pulse using 
Eq. Hwith ^ 1070 cm"\ fin = 1130 cm"\ and 
T = 5 ps. The results of this simulation are shown in 
Figure |3| Figure 13 shows the beat frequency spectrum, 
and original and reconstructured x*^"^-* (fi) Raman spec- 
tra. While the reconstructed lines are broadened, they 
are still quite distinct and would be useful for discerning 
the two molecules. 

NIVI is a flexible tool utilizing ultrafast pulses that 
can measure small or large portions of a Raman spec- 
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FIG. 3: Original and reconstructed Raman spectra in sec- 
ond simulation, (a) Power spectrum of beat frequencies of 
input pulse, (b) Magnitude of Raman susceptibility of hybrid 
DNA/RNA sample, (c) Magnitude of least-squares recon- 
structed Raman susceptibility of hybrid DNA/RNA sample. 



trum of a molecule in a single pulse. It does so by inter- 
ferometrically measuring the anti-Stokes radiation from 
a molecule, stimulated by beats in intensity of an ex- 
citation field. From this anti-Stokes field, the complex 
Raman susceptibility can be estimated. It is especially 
suited to biological imaging because while the pulse en- 
ergy can be large, the peak power can remain small by 
chirping the pulse. For these reasons, we believe that 
NIVI can be a general tool for noninvasively probing the 
molecular content of biological tissues. 
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